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The Kent conductivity x of insulating crystals behaves roughly 
as T * (where T is absolute temperature; with a coefficient which depends 
on the details of enharmonic Interactions and thus is not easy to 
compute. This, project has focussed on the case where K is small indicating 
large anharmonic scattering and correspondingly short phonon mean free paths. 
In this limit the magnitude of x is similar to those found in glasses 
(i.e. 'vl W/raK, to within a factor of 3 or soj . Long ago Kittel^ - suggested 
that the weak dependence of x on T in glasses arises because vibrational 
energy propagates freely (as phonons; over distances (i.e. mean free 
paths , l) not much larger than intermolecular spacings. In this situation, 
anharmonic scattering iG not likely to cause much further degradation of 

and thus i'vconst instead of i'vT \ leaving x'vconst. More recently, 

2 3 

Slack and the author pointed out that even in crystalline material, 
anharmonic scattering could sometimes make i as short as intermolecular 
spacings ,• and that experiments suggest the possibility that in these 
circumstances x behaves much as in glasses. This suggestion is equivalent 
to saying that there is a lower limit, x^^^k^Q/ha below which x cannot 
be driven by any process: alloying, vitrification, radiation dcimage, or 

anharmonic thermal scattering. 

This idea of a lower bound x ^ clearly has adverse implications for the 

ultimate efficiency of thermoelectric power generation. Therefore a principal 

aim of this project was to discover whether experiment supported or denied 

the existence of such a limit. No firm conclusion has been reached. The 

clearest evidence for "saturation" of x at a lower limit comes from 

L 5 

experiments on materials such as CuCl and adamantane . These materials 
quite clearly show x*T 1 at intermediate temperatures but x-*-const'Vi< min at 
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higher T. However, in both cases, the value of k ^ haa a relatively 

strong ahift with pressure. Thie strong sensitivity of to small 

changes in crystal properties suggests that there may still be ways 

available to reduce k , i.e. that is not .an impenetrable lower 

bound. Further evidence is found in a variety of materials of which 

ice^ is a good example, where k in the crystal near T (melting temperature,) 

still has a strong T ^ variation, yet the liquid state value of k is neither 

greatly reduced nor much dependent on T. This behavior seem6 paradoxical; 

K,. should be not larger than k . , yet the solid has k%k . with no sign 
liq min min 

of saturating. Deeper analysis of this situation is inhibited by two 
factors: (i; experimental values of k are particularly unreliable at 

higher T or when k is small, and ( 2 ) theory of k remains very poorly 
developed. Two avenues have been explored with the aim of improving the 
theoretical situation. 

The first avenue is an attempt to provide a simple and reliable way 
of estimating the coefficient of T _1 in the law for good crystals: 
k“A/T. This is explained in detail in the accompanying preprint, ^ which 
proposes a method of estimating T, the'mean scattering rate of phonons 
by anharmonic interactions. From the law r»B/T it was hoped that the 
constant A=const*B could be evaluated. Our analysis shows that A is 
up to an order of magnitude larger than expected from calculated values of 
B. This discrepancy arises from a variety of sources which need a 
detailed anharmonic calculation to sort out. The author plans to do 
such a calculation with R. Shukla next year. 

8 

The second avenue is computer simulation. Mountain and McDonald 
have succeeded in reproducing the law k»A/T by this method in a two 
dimensional case. The author has embarked on such calculations in 
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collaboration with a student, G. Chen, and with D. Emin of Sandia Labs. 

The conclusion so far is that heat conductivity remains incompletely 
understood; further experiments and theories are needed even to clarify 
such a fundamental question as whether a lower limit exists. 
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ABSTRACT 


Vne rate I*, for o phonon to decay by lowest-order anharmonic 
processes is proportional Co T for T>8 D . The Brillouin-zone average 
r of is discussed in detail. An approximate formula for T is found 
which agrees accurately with an explicit calculation for an fee crystal 
with nearest neighbor Lennard-Jones interactions. The Brillouin-zone 
average squared anharmonic matrix element is contained in a parameter 
called which is the cum of the squares of all third order force 
constants. The other Important parameter is the mean square frequency 
which is the trace of the dynamical matrix, or the sum of all diagonal 
(K (£,£)) second-order force constants. The result is a formula 

OCX ^ 

?/ 2 (u z ) 3 where C is (3 tt/ 8) (d^/f^,) and d^ and f^ are 
complicated averages over the harmonic frequency spectrum. To accuracy 
of ^302, C can be replaced by 1. Thus T can be very easily estimated 
when second and third order force constants are known. An "anharmonicity 
temperature" 6^ is defined by the formula 17(^2) “^«T/0^. When T=0^, the 
broadening of the phonons is predicted to be as large on average as the 
frequency, a signal that perturbation theory is no longer valid. The theory 
1 b applied to fee crystals with Lenuard-Jones potentials, and rocksalt- 
structure crystals with Born-Mayer plus Coulomb potentials describing 
anharmonic Interactions, but a shell model describing harmonic properties. 

In all cases 9^ is found to be only one to two times greater than the 
melting temperature. This is compatible with experiment for rare gaB crystals 
but may overestimate the anharmonic strength in rocksalt-structurc ionic 
materials by as much a6 a factor of 2, An average decay rate P 
extracted from experimental thermal conductivity is typically an order 
of magnitude leas than f. 



INTRODUCTION 

» «n ■ 


In crystalline insulators which are not too enharmonic, a thermal 

phonon (u ) dedays ’ primarily through third-order enharmonic coupling 
X 

<V 3 ) , by emitting or absorbing another thermal phonon (u^) . In the 

classical regime (T^O^) the probability is proportional to the thermal 

occupancy of the second phonon, ntw^M^T/fiu^ , giving a decay rate 

increasing linearly with T. This paper explores ways of estimating T , 

defined as the Brillouin-zone average of jj,). The motivation 

1b that T can be estimated in a simple way when information about 

(such as a pair potential) is available, whereas involves difficult 

energy and momentum conservation restrictions. Another motivation is 

that information about phonon decay rates is needed for analysis of other 

processes, especially sound attenuation, and heat conduction. The thermal 

conductivity'*' < Involves a Brillouin zone average 1/T of a reciprocal 

scattering rate weighted by squared group velocities. 

Our search for a simplified formula for F has been guided by analogous 

3 

results in the electron-phonon problem . In particular, the mean electron 
scattering rate, l/x, in a metal with T>9^, is given by ^TiMc^T/fi, where 
the electron-phonon coupling constant X has been extensively studied 

4 

because of its connection with the superconducting transition temperature. 

5 

A formula for estimating X has been developed by Butler et.al, , following 

6 7 & 

pioneering work by McMillan , Hopfield , and Gaspari and Gyorffy : 

X - N(0)<I 2 >/M<oj 2 > (1) 

Our approximation for F Is a close analog of this equation. 


> 

In subsequent sections we derive our approximate formula, test 
it for a nearest-neighbor Lennard-Jonco potential, and apply it to 
rare gas crystals and rocksalt structure crystals. 


2. A, S UM RULE ■ RELATED, TQ,r Q 

We denote phonon quantum numbers by Qi» and more simply, 

O Q 

by 1, The decay rate of a phonon in lowest order is^'’* 

^ *= ti E |V 3 (1,2,3) j 2 {(n 2 +n 3 +l)6Ctim 1 -1ia) 2 -*tia J3 ) 

2,3 

+ 2(n 2 -n 3 )6(tju 1 +tiu) 2 ^Ku 3 )}, 

Momentum conservation restrictions on g 2 and ^ are contained in the 
enharmonic matrix element V 3 (l,2,3). We find it simplifies algebra not 
to exploit translational invariance, but to work instead with general 
harmonic eigenstates denoted by the label i, 

£.,0 

Here i labels the atoms - it summarizes a vector R^° which locates the 
equilibrium site, and an index a or b which denotes the atomic species 
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I 


( 2 ) 

at that site. K^' is the coordinate space dyna ''’al matrix, and 
u fl (i,i) is the normalized eigenvector of the i th mode. The normalization 
and completeness relations are 


£ u a (l,i)u o (£,j) - 6 
£,a J 


(5), 


l o - v»- 


( 6 ) 


The crystal displacement operator 6R Q (£) is given in terms of the 
dimensionless eigenvectors u Q (£,i) by 


6R q (£) » E(V2M £ m i ) 1/2 u a (£,i)A i 


(7) 


*}■ 

where ^ is the dimensionless field operator (a £ +a £ ) and a^ is the 
creation operator. When the states i are chosen to be eigenstates of 
the translation operators, we write u a (*■,!) as 


1/9 ^ *9 * ?£, 

u a (M) - N 1/2 t a (gj.a)e ^ (8) 

where N is the number of unit cells in the crystal and a labels 
the atoms in the unit cell. Whan eq. (8) is used in eqs. 

(5,6), we recover the usual orthogonality and completeness relations for 
the polarization vectors s a . The field operator ^ becomes ^qj =a Qj +a Qj * 

. 'V 'll 'll 

In terms of the eigenvectors u £ , the enharmonic matrix element 


V 3 (1,2,3) is defined by 
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“ < 1/3I > ? ,V 3 (1.2 ( 3) W 3 (« 

V (1,2,3) - Z (tj /2) 3/ 2 (u» w w )“ 1/2 u tt (* » 1) 

3 12 3 0 

x u e (l' f 2)u Y (£«,3)K^(l i i',t») (10) 

“ < K Jl M rV>‘ 1/233E ^ R 2 o 3R £'0 a VY* ' <U) 


A summation convention Is used for repeated Greek subscripts. We are 
Interested in the high T limit of eq. (2) , namely 


r x « (irl^T/fc 3 ) Z |V 3 (l,2 f 3)r| z (tt 1 /w 2 M 3 ) 
2,3 


x I5(u 1 -0) 2 -<D 3 )+26(a) 1 +u 2 -u 3 )] . 


( 12 ) 


Notice that ‘h~ 3 in eq. (12) cancels against "fi 3 in | j 2 (eq* 10); eq.(12) is 
classical. We would like to evaluate F=(3N a ) ^ET^. Instead, we shall 
examine the sum 


A- 5 (8/3N fi 3 ) Z |V_(1,2,3) 
3 a 1,2,3 3 


< V2 W 3 


( 13 ) 


which can be related to f in roughly the same way that McMillan related 
<I 2 > to X. Like <I 2 >, A 3 is surprisingly easy to evaluate. Because the 
factor £D^u 2 aj 3 in (13) cancels against a factor in | | 2 (eq.10), the 
eigenstate labels (1,2,3) appear only on the eigenvectors. The sums on 
1,2,3 are then performed by completeness (eq, 6) . giving the sum rule 


» 
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A - ON)" 1 l (14) 

The quantity »*")] 2 quite generally is short-ranged in | R^-R^ f | and 

JRj-R^Ij even for 1/r potentials, [K^J 2 falls off as r Thus, if 
iB knovm, is easily evaluated, A^ seems to be both a natural and 
a simple measure of anhannonlcity. 

In order to connect A, with f, we define two fairly complicated 
quantities: 


D 


r 

3 


1.2;3l V 3^ 1,2,3 ‘ > 

Z |v„(l,2,3)| 2 w m m, 

1,2,3 J ~ 


(15) 

* 


ii,3l V 3( 1 .2.3)IV^ 3 [MVVV-l-2<(^2- |J 3 )] . ^ 

r l |V 3 (1,2,3) | 2 ti) 1 (aJ 2 u) 3 )“ 1 [6((u 1 -ui ? -a) 3 )+25(m 1 +a) 2 -aJ 3 ) ] 

1,2,3 * “ 

Using these and eq. (11) , we get a rigorous formula for r at high T 


r h (3N J” 1 z r. 

a 1 i 


3wk B TA 3 D 3 /8<(u 1 ^ 


r 


(17) 


The purpose of writing T this way is that A 3 is now fairly simple, and 
the complexities have been displaced into quantities D 3 and <u/*>p which 
we hope to be able to evaluate approximately, by dropping the factor 
| V 3 (l,2,3) | 2 u)^ui 2 u) 3 from numerator and denominator of eqs.(15 t 16). This 
is known as tne "Peierls approximation" (see ref, 1 pp.3S-39). In the present 
context it is somewhat uncontrolled, but will be tested. Then we get approximate 
versions of D 3 and denoted D 3 and 

■ Z 6(ti),-w_-w-)/ Z 1 


1,2,3 


1 2 3 1,2,3 


( 18 ) 


» 
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1.2.3 (i S"V“3 )+24 <“l*V“3 )J 

<(J H > m -J- r U — — — ’ — A 

1 2 3 t l5 ^~ w 2” U 3) +2fi C u i' hl, 2“ W 3^ 


( 19 ) 


F £ On/flJkpT A 3 D 3 /<id* 4 >. 


( 20 ) 


The interpretation of D 3 is that it measures the average decay density of 

states, that is, the number of processes available per unit frequency interval 

for a phonon to decay into two phonons conserving only energy. The quantity 

<w**> provides a measure of the typical value of the factor (l^j) 2 which appears 

in the denominator of eq.Cl6) when the symmetrized numerator usec ** 

■•We expect eqs. (18,19) to be moderately good approximations to eqs. 

♦ « ♦ » • ’ * 

(15,16) not because the weight factor |V 3 (l,2,3j is constant, but 

instead because many states are summed both in the exact forms (15,16) and 
in the approximate forms (18,19). We rely on the cancellation of errors which 
are more random than systematic. In the exact forms (15,3,6) the states are 
restricted by momentum conservation but this is omitted in the approximate 
forms. Of course momentum conservation is very important in eq.(12), and is 
taken into account in the evaluation of A 3 (eq.13) when the exact result (14) 
is used . . . . .. 

It is ‘now convenient to rearrange eq. (20) in order to make several 
dimensionless parameters. First, we introduce the mean square phonon 
frequency 


m2 S (3N fl ) 


-1 


(3N )~ l Z K< 2) (£,fi). 
a £a ua 


( 21 ) 



y 

/ , s 

* 

i 

Hero we ore using oq. (3) and the foct that. Eu* is the sum of the eigen- 
values of Lhc dynnaicol matrix (eq.4) and therefore also equal 

Co the trace of the dynamical matrix. Uaing w* we introduce dimensionless 
versions of the paramntere and < 01 ** > 


S 6?) 1/2 D 3 

(22) 


(23) 


The decay rate is made dimensionless by normalizing to the rms frequency. 

T/® 3/Z " (3«/8) (d 3 /f 4 )A 3 k B T/® 3 . (24) 

3. NEAREST HEIGH BOR LENNARD-JONES MODEL 

To illustrate and test our approximations, we chose a model crystal 
with identical atoms interacting via the Lennard-Jones (LJ) potential 

v(r) - 4c[(o/r) 12 -(o/r) 6 ]. (25) 

For further simplification we let this act only on nearest neighbors (NNLJ 
raocdl) . A rather similar model, but parametrized so as to apply to metallic 
Pb, was treated in ref. 2 by methods somewhat similar to ours. The crystal 
is assumed fee with nearest neighbor distance chosen to minimize v(r) 
in eq. (25), i.e. where v(r Q ) — e. To evaluate 


J7 and IV_ we need expressions for the derivatives of the additive energy 

J i 

of pair potentials: 


&2p (R« rt "^o 

b i a " t,,5 °e 


(26a) 


a^E , ( R t.n~ R t'f^ ^ R tB~ R t'8^ ^ K tv~ R t'v^ 

Vve 3 h'y " ° 3 ^ l R i“ R jf 1 3 

*V *Vi 

. , t ... (R ia _R £ , a ){ 6v +(R t6" R t'6 )5 vo +<R tv' R t , Y ,i aS 

+ a 2<« > IvvT 

fV 'V 


(26b) 


These expressions are valid for V+l. When the corresponding expressions 

are: 


' 2? 32E/SR «a 8R 1 '6 <27a) 

33e/3 W R *y ‘ 33E/3R ta 3R H'6 SR J'Y- <27b) 


The coefficients b^,b 2 i a 2 * a 3 are 


, 1 5v 

D 1 r 3r 

, 3 u 

b 2 ’ 1 3? b l 

a 2 (r) - |j b 1 

a 3 (r) - r> |j: <± ^(D). 


(28a) 

(28b) 

(28c) 


(28d) 
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From these formulas one can get a general expression for I? (eq, 21) and 
(eq. 14) for any material described by pair potentials (PP) 

I?(PP) - ON) -1 Z M^Ib-U-Jl'HSb. (£-£')] (29a) 

tyl' * 2 *■ 

A 3 (PP)«(3N)"^, (M Jl M^)" 1 [a 3 (A-£’) 2 +6a 3 (£-t')a 2 (t-ll , )+15a 2 (£-t , ) 2 ]. 

(29b) 

For the NNLJ model these become: 

I?(NNU) » (4/M) [b 2 +3b 1 3 » 288c/2 ly ' 3 Mo 2 (30) 

A 3 (NNLJ) “ (12/M 3 )[a 2 +6a 3 a 2 +15a 2 J - 13 ,903,488e 2 /M 2 o 6 (31) 

Equation (24) then becomes 

Tf (to?) 1/2 (NNLJ) - (3a/8)(d 3 /f 4 ) (149/128) (k B T/c). (32) 

The parameters d 3 and were calculated to 1 % accuracy by numerical 
evaluation using a tetrahedron program. The values obtained are d 3 =,095, 
f 4 “,129, and .737 . As a test, 1? was found to be 287.0(e/2^ 3 Mo 2 ) , 

agreeing well with the exact value of eq. (30). The corresponding value 
of T/6?) 1/Z in the NNLJ model is 1.01(k_T/e). 

a 


Our estimate of f involves the uncontrolled approximation of replacing 

i 

| V 3 1 2 u 'i u, 2 a> 3 ky 1 * n going. from (15) and (16) to (18) and (19). To test 

this we have evaluated F directly from eq. (11) using the correct frequencies, 

polarization vectors, and matrix elements of the NNLJ model. Crystal 

momentum conservation was explicitly included, but the energy-conserving 

6 function was replaced by a Lorentzian Im(x-i6) ^/n of width comparable to the 

finite mesh size increment Aw»|dw/dQ|AQ, Using 4000 k-points in the 

- — 1/2 

Brillouin zone sums, the answer was r/( uJ 2 ) «l,08(kgT/e) . This anewer 

was stable to about 5% under changes in mesh size and 6 . The good agreement 
with our approximate answer, 1.01(kgT/c), exceeds reasonable expectations 
and must be fortuitous. We do not expect the accuracy of eq. (24) to be 


better than 20-30%. 


4, FULL LENNAR!)- JONES MODEL 


The purpose of truncating the (6-12) potential at nearest neighbors 
in the previous section was only to reduce the computer time needed to 
calculate the matrix element in the exact calculation. Our approximate 
formulas are as easy to' evaluate with gll neighbors as they are for first 
neighbors only. We need the lattice sums 

Z £ E d n /|i| n (33) 

n * 

where d is the nearest neighbor distance and i runs over lattice vectors. 

For an fee lattice, the values of Z needed here are Z o b 12.8019, Z_ .«1?.0590, 

n o id , 

Z l8 «12.0130, Z ^-lZ, 0015, and Z 3Q -12.002. Then and w 7 can be evaluated 
from eqs. (30,31): 

T? “ (8c/Mo 2 )p 8 [22p 6 Z 1 4-5Z 8 ] (34) 

eA 3 - (8c/Mo 2 ) 3 (18p 18 )[8575p 12 Z 30 -2716p 6 Z 24 +220Z l8 ] (35) 

where p-o/d. In classical approximation at T-0 the atoms are stationary 
and minimize the total energy. This occurs at p-p - 0 / d where d /o«(2Z. 0 /Z,) 
-1.0902. At this value of p, the dimensionless anharmonic parameter 
eA 3 /(o> 2) equals 0.7172. The values of (d^f^) have been evaluated to 1 % 
accuracy with a tetrahedron program. Results are shown in column 1 of 
table 1. The widch-to'-f requency ratio r/(w 2 )’ L ^ 2 is 0. 76(kjjT/e) , 25% 
smaller than in the NNLJ model. This is still a remarkably large number. 

Rare gas crystals have melting temperatures T^O. 7 (c/k^) at 1 atmosphere. 

In our notation r«-2Im£ is the full width at half maximum for a Lorentzian 
lineshape. Thus the average phonon width at T=T m is predicted to be ^50% 


of the rms frequency. This is so large that 1 , at least for a significant 

subset of the phonons , 2nd order perturbation thdory can no longer be 

accurate and quite probably we are outside the radius of convergence of 

the perturbation series. Bohlln^ found by direct evaluation of eq. (2) 

. that LA phonons in Ne at T-4.7K (about 25 % of T ) had widths T as large 

as 40-50% of the frequency. The failure of second-order th '•.rmodynamic 

perturbation theory for T^T^/3 had been noted by Klein et.al.^, and has 

12 

recently been examined to higher order by Shukla and Cowley. Neutron 


13 

experiments in Kr have seen values of comparable -to for zone 
boundary LA phonons at T' close to T m . Molecular dynamics simulations 


14 


of S(Q,w) for LJ systems have also seen broad zone boundary LA response 

functions near T m . One can then ask whether second order perturbation 

theory gives qualitatively correct trends even in the regime where 

the justifications for perturbation theory fail. We are not able to answer 

this quantitatively, but published dispersion curves, lineshapes , and 

simulations all suggest that ill-defined phonons with Pvuj are rarer than 

our estimate gives. In other words, the actual behavior of the strongly 

anharmonic system tends to give quasi-particle-like response even when 

perturbation theory says that the quasi-particle picture should no longer 

be valid. The other possibility is that ill-defined lineshapes are less 

likely to appear in publication than well-defined ones, and that our 

estimates remain reasonably accurate even near T^. 

The Bource of the large anharmonicity lies in the steep and one-sided 
-12 

nature of the r potential used to model the large repulsions when closed 

shells overlap, combined with the softness of the potential for r ^ r m ^ n * 

15 

These factors also cause a large thermal expansion of ^3% at T ffl which , 


(significantly alters the phonon response at higher T, In fig. 1, i s 
plotted versus d/a, showing a dramatic downward shift when d/o increases 
by 3%. Thus It is important to use the corrected harmonic frequencies 
at temperature T (quasi-harmonic model), and the corrected anharmonic 
.matrix elements. The measure of anharraoniclty also decreases drama- 
tically as d/o increases, but not as rapidly as (u7) 3 decreases, so that 
the dimensionless factor cA,/(m^) 3 is quite strongly increasing as d/o 
increases, as seen in fig. 1, ’-We have recalculated all parameters 

at d/o«1.12, and the results are in column 2 of table 1. The width-to- 
— 1/2 

frequency ratio r/( u 2) becomes 1.26(kgT/e), 66% higher than at 

d/o»1.09. Thus using the quasi-harmonic approximation as a basis for 
doing perturbation theory only makes the anomalous magnitude of 
17 (w 2 ) 1 ^ 2 more serious.. 

Finally we turn to thermal conductivity k . From Boltzmann theory 
we obtain 

* " < VT >' 1 ^,HVQQ'“Q' v Q' J . < ‘ Sn Q /a “Q' ) (36) 

where V is the volume, Q is short for phonon wavenumber and branch Q j , v 

f\, Qx 

is the group velocity 3io / 3Q , and n is the equilibrium Bose-Einstein 

y x 4 

distribution. In relaxation-time approximation, the scattering operator 
(k) 

Sqq, is rj 6 qqi where the superscript tc reminds us that this differs 
somewhat from the quasi-particle scattering rate especially in that 
N processes (non-Uraklapp) are not fully effective in damping the heat flow. 
At T>0jj, n^ is k B T/ti(UQ and (36) becomes 
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Thus we define 


V ‘ (3N / l jjb v 8 J 

(38) 

l/r W £ Itv/Nkj V 7 

(39) 


where N is the number of atoms. The experimental value of K for Ar at 
high T can be expressed as"^ 

1/k - (5.3xlO" 2 )M 1 ^ 2 d 2 c" 3/2 T. (40) 

-2 

The number 5.3x10 comes from constant volume measurements by Clayton and 
Batchelder.' 1 ' 2 No other rare gas solids have been measured systematically 
at constant volume, but the trends suggest that eq, (40) should be reasonably 
accurate for all of them. 

We have calculated v 2 “24.5e/M for LJ crystals at d/cf“1.09 and v^«l4.2c/M 
at d/o»l,12. Using the former value as more representative of the constant 
volume conditions of ref. 16, and the value oi^ , =449.2 e/Md 2 obtained from 
eq. (34) at d/o«1.09, we find 


r^/Otf 2 )^ 2 - .087 (kgT/e) . 


(41) 


(k) — 

Thus r w is less by a factor of 9 than the theoretical value of T. 

Three causes contribute to this discrepancy. (1) Since T is 

defined by an average of l/T^, it weights small values of most strongly 

(occurring for small Q acoustic phonons). This is reinforced by the 

weighting factor of v* in eq. (37). (2) Umklapps contribute fully to 

Q 


(ic) fiO 

Tq ' and Tq, but N-processes occur more weakly in 


This should 

contribute less than a factor of 2 to the discrepancy, (3) As previously. 


17 


. I • ‘ 

« 

mentioned, quasiparticles may be better elementary excitations than 

perturbation theory says they should be. This idea is supported by the 

data of ref. 17. In cases where quasi-particle approximation is known 

18 19 

to fail, k seems to saturate * at a value K min rather than decreasing 
as T _1 . The date of ref. 17 obey k«T - 1 quite well. 

It is not possible without detailed calculations to further subdivide 

( K ) *** 

the cause of the large difference between V and T. The most detailed 

20 

calculation to date seems to agree well with ref. 17 but not to shed 


much light on this question. 
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5. UOHN-HAYEK MODEL FOR ROCKSALT STRUCTURE 


Rockaalt (Nad) structure compounds hnve been the subject of much 

theoretical work, and there exist data on phonon lincvidths In NaCl 2 ^" 2 ^ 

‘ 27 

and on the thermal conductivity of several compounds in the group. 

Anharmonicity in these materials at T£600K is weaker than in the rare gas 

crystals near their melting points, allowing greater confidence in pertur- 

28 

batlon theory. Our starting point is a Born-Mayer-type pair potential , 

consisting of a Coulomb term and a nearest-neighbor repulsive exponential. 

It is well known that the phonon dispersion is not very well fitted by 

such potentials, especially the optic phonons, but it is believed 2 ^ that the 

anharmonic part of the interatomic force is adequately treated in this way. 

Thus, we will calculate from this pair potential, and take quantities 

like the mean square frequency and group velocity from shell models that have 

29 

been previously fitted to detailed spectral data . This model potential 

has the advantage of permitting us to work out closed form expressions for 

most of the interesting quantities, and to make direct comparison to other 
2k 

calculations . The Born-Mayer parameters will be taken from standard fits 
to the lattice constant and compressibility; one could trivially extend our 
results to a three-parameter model by introducing a non-integer effective 
charge, as when, for example, one also wishes to fit to the total binding 
energy. 

For a rocksalt-structure crystal in which the atoms carry charges ±Ze, 
the pair interaction is taken to be 


z * z r - g2 

r w 





“ r U’ /p 


+ C e 


(42) 



I 


I 

• 

The second term approximates the overlap repulsion between adjacent 

atoms, with C, p being chosen to fit a given compound. This term is 

assumed (as part of the model and not as an additional approximation) to 
* 

b 1 ' nonzero only for nearest neighbors. The cohesive energy per particle is 
given by 


U " zT S (43) 

£, A* 

where the sum is over all sites of the lattice except that Now let r 

. o 

be the nearest neighbor distance, equal to a/2 where a is the lattice constant. 
Cutting off the second terra at nearest neighbors and introducing the 
Madelung constant . o»l. 744 ... , we have 


U(r o ) 


aZ 2 e 2 . o 

+ 6C e 


r Jp 


(44) 


The lattice is stable for that value of r which satisfies 

o 


aZ 2 e 2 r - * ~ r - /p 
6Cp 'p 


f 0 \ p o 
(— ) e 


(45) 


Taking derivatives of $, we calculate the quantities of (28a, b): 


1 1 3e 2 Z.Z„, r i 1 j i /p 

a (£,£.») o — — rfi" ± — . — ^ At 


r £Jt ’ 




n 


(— — + — )e 
* pr JU' r U* P 


a {*, £')“♦ m - — — 4>" + - £t — ■ - -(^ + -5 — + — ) e 

3 r n' nv r u’ *u' pr n' 


(46a) 
it' 


_ ^(1 B » /P 


(46b) 


where it Is again understood that when lattice sums are taken, the 

exponential terms only include nearest neighbors, We can also work out 

a Born-Mayer expression for the mean square frequency using eqs. (28,29) 
• \ 


IF 


(i- + i_)£(i _ 2-) r " r ° /P 

'M, M Vp rj 


le ' . (47) 

•+ K »* ‘o 

By the relation (45) all long-range terms are eliminated from (47) , which 

.30 


becomes a purely nearest-neighbor quantity*, this sum rule" 
of the coulomb potential's satisfying Laplace's equation. 


is simply a result 


We next evaluate Ag, starting from expression (29). The lattice sums 
over the long-range Coulomb terms need to be done with some care; the 
result is 


a - 3 ,1 . l» q5Z**e 4 _sc l2Z 2 e 2 C 

3 " M.M 'M. + M 'I. r 8 a ~ r ** P 

+ - + - r o l 0 

,1 . 3 + 3 . “ r o /p C* 1 6 ,_12 „ . 

( ? r+ r n p + rS )€ 


6 ,/V^ 


45Z t4 e* 4 Z B 


fee 


16 


(1— i_ ) 2 ( JL + i- ) 

V M ' V M, M ' 

• 4 - + 


(48) 


where Z 


fee 
8 


SC 

■12,8019... and Z^ ■=6.9458... are the lattice sums defined in 


eq. (33), but specifically for face-centered and simple cubic, respectively. 

g c 

Since the last term turns out to be small (for realistic parameters) and Zg 
is only 16% greater chan the nearest neighbor value, 6, we see that is 
dominated by the nearest neighbor force even in the case of long-ranged Coulomb 
potentials . 

It is instructive to compare our with the same quantity calculated from 

24 

detailed models that closely fit the experimental . Using the shell model 


for NaCl, for example, we have obtained the density of states F(w) and 

from it calculated The result is tow z ^ Z /kg“212.5K, while our 

Born-Mayer model gives li^ 1/,Z /k n -225.6K. This 6% discrepancy in the 

rtns frequency becomes a 40% discrepancy in <uj 2 ) 3 which is needed in eq. 24, 

The agreement is less good, if we examine fu* 2 , which characterises 

the shape of the spectrum. The NaCl Born-Mayer potential gives 0,291, 

while the numerical result from the shell model is 0,523. The large error 

is not surprising in view of the well-known underestimation of optic 

frequencies in the simple model, Thus, it is preferable to use F(w) from 

the full (shell model plus tetrahedron program) calculation to get the 

factors d 3 , f^, and ^ i n (24). It is no more work to do this, since 

even with the Born-Mayer potential the full F(uj) (as distinct from its 

first few moments) would have to be calculated numerically to get d 3 and 

f^. In the process of doing this we have also numerically calculated 

reliable values for the mean square phonon velocity v^" of (38). 

We have obtained results for LiF, NaF, NaCl, KC1, KBr and MgO, these 

being representative (highly ionic) rocksalt-structure materials for which 

good thermal conductivity measurements have been done; phonon linewidth 

21-26 

measurements appear to be available only for NaCl. In Table II our 

value of for NsCl is compared with the measured ratios for 

several phonons. Our results for to 2 , A 3 , v^, and f are given in Table III, 

together with the experimental thermal conductivities of all six materials, 

(O 

and the values of T' ' derived from them. 

For NaCl, our value of is bigger by factors between 1.2 and 7 

than the various experimental ratios. Unlike the case of argon, one cannot 

ascribe this disagreement to a breakdown in perturbation theory, 6lnce our 

predicted relaxation times at room temperature are an order of magnitude smalle 

24 

than tHe corresponding frequencies. Eldridge and Stahl , with a similar pair 


n 


i t • 

potential ami shell model, obtained phonon liiicohnpea in reasonable agreement 

with experiment. Thuo we would expect our value of r/(uJ7)^ 2 to be rcliab’ 

» 

i 

It io conceivable that many phonono, an yet unmeasured, may have large values 

of Tq. Another possible explanation is that our neglect of momentum conservation 

is a poor approximation to make in the decay of the optic phonons, which have large 

regions of flat dispersion where energy conservation is easily satisfied , 

(k) 

Table III also shows scattering rates I" derived from measured thermal 
conductivities x using eq, (39). The values of F s are almost all an order 
of magnitude or more smaller than the theoretical 7 values, similar to the 
case of rare gas crystals. Six possible causes of the discrepancy ore; 

(1) possible inadequacy of the Born-Mnyer model; (2) inaccuracy of the 

approximations leading to eq. (24 ) ; (3) failure of perturbation theory! 

’ — 

(A) genuine differences between P and T arising from the suppression of 

Uraklapp scattering 32, in ; (5) genuine differences between and T 

arising from the bias in eq. (37) toward long-wavelength acoustic branches 

(r) 

with large Vq and small 1^ > (6) experimental vmcertainty in x, especially 

from possible failure to subtract radiative transport. 

We believe the differences are genuine. Umklapp's (cause no. (A)) 

probably account for a factor of two and most of the rest is cause no. (5) 

- the variations of P. with branch Qj are quite extreme and different methods 

Qj n. 

of averaging can generate an order of magnitude difference. This can be seen 
by a study of table II and comparison w'< th F / (m*) ^ in table III. Reasons 
(1-3) are ruled out by the success 23, 2 ^ of anharmonic perturbation theory 
based on Born-Mayer potentials for individual widths Tq j , and by the t^st 
we performed on our approximations in sec. III. Experimental accuracy is 
alwavs a problem in measurements of x, but this is likely to play only a 
mitsvr role here. 
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6. CONCLUSION 

It Is appropriate to emphasize the virtue of P/(t?)^ 2 as a measure 

of anharmonicity. Unlike the usual measures, i.c, Criinciscn parameters and 

* 

thermal expansion, which measure long-wavelength anharmonic effects, P probes 
all wavelengths democratically. The fact that T docs not accurately predict 
any particular width P^ and that P overestimates the heat conduction scattering 
rate P^ does not invalidate our argument. The ability to estimate T easily 
may provide a new perspective on the problem of anharmonicity , and should 
accurately indicate the adequacy or inadequacy of anharmonic perturbation 
theory. 

A convenient way to reexpress the information in P is to define an 
'’anharmonicity temperature" 0^ as the temperature where T is as large as 

-r 1/2 

(wO and perturbation theory falls: 

r7(w*) 1/2 s T/e A . (49) 

This definition assumes that T>0^ and that perturbation theory in lowest 
order gives the dominant behavior. For the Lennard-Jones crystal, 0^ is 
1.32e/k B based on the zero temperature nearest neighbor distance, 1.09o, or 
0,79e/kg based on a high T distance, 1.12c. Thus is higher than the melting 
temperature T M ^0.7e/kg, but only by a factor *1.1-1. 9. Similarly for NaCI 
structure, tnble 1X1 shows that 0^ is typically 3-2 times larger than T^. 

The highly anharmonic nature of these materials when T/feT^ is not widely 


appreciated. 


r 


o 


Finally, given the 20-30IS uncertainty we ascribe to our eq. (24 ) for 
r/(w 2 ) • the difficulty of evaluating and f^, and the unpredictability 

of Fqj and it becomes appropriate to offer a simpler formula. The 

factor d 3 /f^ is given for various cases in tables I and III, From these 
numbers we can expect d^/f ^ to be moderately insensitive, to details. The 

factor (3n/8) (d 3 /f / , ) in eq. (24) can be replaced by 1 with an error typically 

* 

±20£; in the extreme case of KSr, the error is a factor of 2. Then eq. (24) 
is replaced by 


r/U*) 172 = T/0 a a A 3 k fl T/ <0,2)6, (50) 

The parameters of this formula, and w 2 *, are numbers which can be estimated 
on the back of an evnelope when a model is available, and provide a 
surprisingly simple and accurate way of characterizing anharmonicity . 
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TABLE I . Calculations for Lennard-Jones Potentials 
of the Parameters Entering Eq. (24), 


Full LJ , NN Model 


d/o 

1.0902 

1.12 

1.1225«2 1/6 

eA ^/(w^) 3 

0.7172 

1.2774 

1.1641 

d 3 

0.127 

0.114 

0.095 


0.141 

0.136 

0.129 

V f 4 

0.90 

0.84 

0.74 

II 

r IG ?) V1 

O.76k 0 T/e 

1 , 26k_T/c 

D 

1.01k o T/e 

D 



28 


TABLE II . Linewidths of Measured Phonons in NaCl 
at 30QK. Our Theoretical Value of 

at This Temperature is 0.186. 


Phonon 

q 

LO 

0 

LA 

Is- C 1 - 0 - 0 ^ 

LO 

several 

TO 

0 

TO 

0 

TO 

0 


*at 290K 


P /o> References 

_g a 


o,0,ll 

21 

o-0. 05 

22 

0-0.16 

23 

* 


0.025 

24 

0.04 

25 

o-0. 04 

26 



TABLE III ; The first 5 rows are calculations from shell-model phonon parameters (ref. 29) of 
rms frequency, group' velocity, etc. Row 6 gives calculated from Born-Mayer 

potentials (eq.48). Rows 9jl0 compare from theory and from thermal conductivity 

k. Experimantl values of tc (row 11) are from ref. 27. 
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(W/mK) 20. 19. 6.2 7.0 4.7 *u60 

(expt, 300K) 



FIGURE CAPTION . 

Fig. 1: The lower curve is the hennard-Jones (6-12) potential versus 

r/o for a range of separations near the minimum at r /a m 2^^ 

tain 

i 

*1.1225. The upper curves give dimensionless measures of the 
mean square frequency and the anharmonic parameter A 7 /(w?) 3 

w 

* versus d/a where d is the nearest neighbor spacing. In 

classical approximation at T»0, d/a takes the value d 0 /o»(2Z 

° 12 6 


1.0902. 


* f> +r , 
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